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SUBJECTS PROGRESS REPORT, AUGUST 9 THROUGH SEPTEMBER 5» 195fr 

Tos Jay ¥. Forrester 

From: Scientific and Engineering Computation Group 

1. MATHEMATICS, CGPIJG AMD APPLICATIONS 

1*1 Introduction 

During the four week period covered by this report, 722 coded 
programs were ran on the time allocated to the Scientific and Engineering 
Computation (S&E.G) Group represent-in j work that has been carried on .in '41 
problems ♦ 

Section 1.2 contains progress -.'■..-/ports && submitted by the programmers 
together with an indication of the :■;;■ hiine tiiis (and, in some cases,' of the' 
programming time) expended on each ;;robleia. 

Two simulated computers were developed by the S&EC staff and used by. 
the 65 students enrolled in the si:...* ..: &u.:.<.r Session Course 6.531 on 
^Digital Computers; Business implications® This course was given under 
the supervision of Prof, C. M c Adaru:; and Dr S e Gill during the period 
August 16-27. -he problems that were solved by the students on these 
simulated computers are described briefly under S&EQ proble:.: 19B. The 
"computers" themselves are described under problems 196 ana 197. 

Detailed initial descriptions are provided below for problems 195 
(analytical study of intestinal motility records) and 207 (solution of 
loth and 19th order differential equations for the M.I.T. Flight Control Lab), 

During this period the S&EG Group played host to 220 visitors. A brief 
summary of the groups represented by these people is given under problem 131 T 

1.2 Programs and Computer Operation 

100 a .ThgLi?^ Demur jian, 25 hours; 

Frankovich, 16 hours; ' Helwig* $0 hours?; Siegel, 4 hours; Xfetson, 
56 hours; Porter , 2 hours; W3I, 719 minutes 

CS II has been modified to operate under the guidance of ra director". 
control tapes , These are Flexo tapes containing information equivalent to 
that in a performance request and typed in a manner to be described in a memo 
being written hy F. Helwig. These tapes are read by .the new Input Program 
(soon to be recorded on drum group 11) via the mechanical tape reader. 

J.'-K. Frank ovich 
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A director tape for a particular run is a paper tape which is obtained 
by typing the performance request for the nan on a Flexowriter. Under certain 
restrictions it will soon be possible to place the director tape in the mech- 
anical tape reader 9 the tapes involved in the run (spliced together in the 
proper sequence) in the PETR, and then to execute the run by a single pushing 
of the read-in button. 

A director tape program which reads in words from the director tape and 
effects the appropriate "button pushing" has been written and checked out. 

The introduction of director tapes requires that several of the S&EG 
recorded utility programs be modified. Some of these mods have been written 
and tested. Work continues on those remaining. 

A memorandum describing the standard abbreviations to be used in filling 
out S&EC performance requests is being written. This will also contain a 
description of the director tapes . 

F. Helwig 

Two (30 9 0) divisions for quantities in fractional form have been complet- - 
ed and are being written up. 

One routine using the binomial iheroem and 84 registers produces results 
with a minimum accuracy of 2? digits . The other routine using 140 registers and 
the method used in the Summer Session integer divide routine gives a result 
accurate to the l$st digit. 

R. Mats on 

As mentioned in the previous progress report the routines for use in 
automatic curve plotting on the scope have been recorded on magnetic tape. 
We are now in the process of testing the conversion of tapes using these 
instructions and of testing the routines themselves, 

I.- Saber 

106 C. MIT Seism ic Projects Grine, 50 hours 1 Simpson, 50 hours; WM, 

466 minutes 

The time used during the past 4 weeks, by problem 106 represents 
almost exclusively experiments conducted on multi-trace seismograms not 
utilising the linear operator as such. Instead, various new mixing tech- 
niques are being tested to derive information supplementary to filtering. 
In one experiment Fourier Transform phase angle mixing is being used to 
determine input signal-to-noise ratio as a function of frequency.. Both 
the theoretical and practical aspects of this problem do not permit analytic 
solution and are being computed. One example tested under controlled con- 
ditions yielded expected results, but full interpretation must await the 
theoretical solution. 

A correlation function mixing technique being tested should theoretic- 
ally yield time varying estimates of signal power with dip of reflecting 
horizons as a parameter. Preliminary tests have yielded significant in- 
formation from a set of records previously considered worthless. 
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120 Do The Aerothermopressor s Porter, I hours Gavrils 101 1 124 minutes 

A study of the fourth-order Runge-Kutta method for solving the 
differential equations of the aerothermopressor has continued and is nearly 
completed. This study involved carrying out the computations for a typical 
case using the Runge-Kutta method with various increment sizes and compar- 
ing the results with those obtained from the Euler method for various increment 
sizes , 

It has been found that the spurious behavior in the Euler solution which 
had been believed to be the result of truncation error (see progress report 
of July 12, 195*0 has been completely eliminated by the more refined numerical 
procedure. For increment siaes of 1$ (i,©. the range of the independent 
variable divided into 100 parts), the final computed pressure was in error 
using the Runge-Kutta method by an amount equivalent to that which would arise 
with the Euler method using an increment of .0^$, For this accuracy (within 
about .3/6 in the pressure) the Runge-Kutta . method is therefore about six 
times as fast as the Euler method. An increase of increment size to 2$ win 
the Runge-Kutta method was found to give acceptable results with the- error 
only twice as large as that corresponding to the increment size of 1$. 
The apparent linearity of the error using the fourth-order Runge-Kutta method 
is not yet clearly understood, and further computations are being carried 
out . 

The results so far obtained s&em very promising from the point of view 
of reducing the time of computation by using large increments over part of 
the range of the independent variable. It appears that the computations 
must begin and end with 1$ increments s while in the mid-range, steps of k% 
may well be taken* 

Previous plans to write a completely new program remain unchanged, aid 
this work will begin as soon as the necessary preliminary details are worked 
out in the present study» It is planned to use the current subroutines to 
continue. aerothermopressor computations to a limited extent while the new 
program is being written. 

126 Co Data ReductionsD .T.Ross ;¥,¥olf ;D. Hamilton i Ml, 153 minutes 

' This progress report covers the last two monthly periods since vacation 
schedules interfered with the writing of reports, Three complete versions 
of the Data Reaction Program have been tested with and without the Mistake" 
Diagnosis Routine (MDR), described previously, which was developed under 
this problem for this purpose. The MDR has served its purpose very well 
and its use is now a routine matter. The first version of the Data Reduction 
Program was a simplified version and has now been discarded since it has 
served its purpose of simplifying the debugging process. The second version 
of the Data Reduction program has been completely checked and has given 
useful results The second version will also be discarded, however, in 
favor of the third which uses modified computing techniques for higher 
speed and greater efficiency . Checking of the third, the so-called Basic 
version, is progressing rapidly with the use of the MDR. Further modifications 
of the Basic version will be made to compute different answers. Each successive 
version will be an extension of a previous version, designed to handle more 
complicated problems. 
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Work is progressing on another phase of this problem , using auxiliary 
in-out equipment both as an integral part of the Da^a Reduction Programs 
for monitoring purposes and as an aid to experimenting with improved 
' "computational techniques . Intervention registers and scopes -will constitute 
the' ''working parts of this system with efforts being made to provide the 
flexibility, normally associated only with analog computing programs. Present 
work is concerned with the development of programs to decode intervention 
register inputs and an alpha-numeric scope routine equipped to handle (24-, 6), ■ 
(I5t>0) and alphabetical characters with decimal points suitably placed, etc. 
A large-letter title display is also being written* ! 

131 . Special Pr oblems (Staff Tra i ning B Demonstrations „ etc.) ; Demur jian, 
6 hours 

During the past period 220 people attended the tours of 'Whirlwind. 
On August 12 the regular monthly demonstration of the computer facility was 
heldo There were 33 visitors present .' 

Sixty-eight students from the Summer Session course 6.531 Digital 
Computers? Business Applications were shown the computer and control area 
on August 18 and 19 in preparation for their later use in solving assigned 
problems . 

On August 23 the movie ' 3 Making Electrons Gount 55 was shown to the teachers 
attending the M. I. T« Mathematical Institute Meetings. Sixty-eight of these 
teachers visited the Digital Computer Lab. for a demonstration on August 24. 

The last demonstration of the period was for 51 students from the 
Servomechanisms Laboratory at M. I. T. on September 1„ 1954. 

132 C-. Revisi on ju-Exteng [ion ,_and. Testing, of the Subrou tine Library Used in Programs 

for Obtaining Data for the Humerica J-ly Controlled Milling Machine s 
Runyon, 20 hours? ~%MI 36 minutes 

Time used this month was devoted entirely to milling machine tape 
checking o A program modification for printing milling machine order totals 
at preset intervals .was used to facilitate location of errors.. 

141. S&EC Subroutine Study 5 . Mahoney, 49 hours? WI, 80 minutes 

Several modifications or replacements of existing subroutines have been 
made and are now in the Library of Subroutines. These are 
LSR Ful _ Exponentials subroutine 
Fu2a Square root (24 9 6) 

Fu2b Square root (30-j 9 j) (four iterations) 
Fu5 Hyperbolic functions 
Mkk- Matrix diagonalization 

The Auxiliary Buffer Routine (LSR SPl) has been obsoleted since it is un- 
necessary with CS IX in operation 
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143 Do Vibrational Frequency Spectru m of a Copper C rystals Jacobsenj 
WHI, 97 minutes 

The previous program developed in problem 1^3 calculated the vibrational 
spectrum for a face-centered cubic crystal. This program is now being modified 
in order to calculate the vibrational spectrum for a body-centered cubic crystal 
and will include 1st, 2nd, and 3 r d neighbor effects as did the original program. 
In particular;, the calculation will be made for c<- iron employing the tensor 
force constants determined by H. Curien, College de France , 1952. This calcu- 
lation will conclude problem 1^3 « 

' '' 

156 Ao Evaluation of the Reflection Coefficient in a Semi°?Infinite Rectangular 
Wave Guides A, Balserj WMX, 8 minutes 

Another test run was held for the range Q^^^TT The results 9 
when compared with check results which are available for some values. of &< , 
were off in the second or third place . Consequently , a test run was then 
held f or <=*- ~ and a post-mortem was made. The results are now being studied 
in an attempt to detect the error* 

159 D Water Us e in a Hydroelectric Systems Littles 1ML 9 107^ minutes 

Improvements in the optimization program and a graphical scope output 
routine are now incorporated into the main tapes. Final results are being run 
off for the model of the hydro system under study. The scope output is excellent 
for seeing how the system, using the operation determined by the program, be- 
haves when the historical river flows are the input. 

166 Co' Const ruction an d Testing of 3eltajirig>: Flutter Model. Chen, 3 hours J 

WWI, 95 minutes 

During the past period, several runs were made to find the Ac's'. 
It is intended to make a parametric study in order to determine the effect' 
on £s>e s s due to a change of x. The existxng tapes can be used for this 
purpose. 

167 Bo Transie nt Effec ts in, Distillations. Jordan, $0 hours? O'Do^nell, 

160 hours 5 Smith, 5 hours s' Porter, 2 hours? WI, 1013 minutes 

About half the machine time was utilized in obtaining additional 
solutions for the case of a transient in a continuous distillation column 
caused by a step change in feed composition. Two correlations had been 
obtained for the data available previously. The additional solutions were 
used to determine the effectiveness of the correlations in predicting results 
for different parameter combinations. The work showed that the correlations. 
were quite successful for conditions within the range previously studied. 
However 9 some extreme sets of parameters, outside the isnge studied, gave 
results that were not predicted by the correlation satisfactorily for engineer-, 
ing purposes. This was rather probable and further work will be done to improve 
the situation. Jordan has completed his Master's thesis on this subject and 
has submitted it in the Chemical Engineering Department. A considerable im- 
provement in understanding of the situation has been achieved. 
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The remaining machine time has been spent on the problem of batch 
distillation take-off. Most of this time has been utilized for obtaining 
solutions for the situation described in the Progress Report of August 8, 195^. 
This is the limiting case of batch distillation take-off assuming the hold-up 
is negligibly small. The basic relationship is s 

V - exp 

where ¥ is the fraction of the charge not yet distilled x^ is the distillatp 
compositions and x r is the still pot composition. From a combination of mater- 
ial balance and equilibrium relationships x can be expressed explicitly in 
terms of x~. However the only way x n can be obtained irom x is by iteration. 
The initial value of x-. corresponding to x is calculated by iteration. There- 
after Xq is chosen and x calculated/ The resulting increments in x_ are un- 
evenly spaced'. One of the problems was the question of how to carry out the 
integration, 
t 
Initially it was hoped that using the trapezoidal rule 8 solutions of 
sufficient accuracy could be obtained with reasonable machine time, .By run- " 
ning two sample cases on the ma chine s each with at least three different interval 
sizes 9 it_was determined that the results 9 «" 6 for any value of x when plotted 
against h where h ? is the increment size 9 was a straight line. This implies 
that Richardson 9 s h extrapolation might be used. Briefly the scheme actually 
followed in the solution was to calculate values of l f J using two different 
interval sizes alternately with one half the other. That is 9 two steps would 
be taken using the smaller interval size and then one with the larger 9 at the 
end of which there was available "H at an x Tt calculated using two interval sizes. 
Then the mathematical equivalent of Richardson's h^ extrapolation would be done 
to get a corrected value. This was printed out. The results were very satis- 
factory.. Less machine time and output was used than would have been necessary 
if the two calculations had been done separately (because it was not necessary 
to calculate x from x^ again f>r the larger step). Also less machine time was 
used than would have been necessary if each calculation had been done only with 
one increment size small enough to insure sufficient accuracy, A check was 
made at the time the extrapolation was carried out to see that it was not too 
large. 

Some machine time was spent getting additional data for the case- of 

batch-distillation take-off with holdup. 

At present, most of the effort is going into attempts to understand and 
correlate the data for batch-distillation take-off, Scr.ie machine time was 
used attempting to explore a correlation which involved solution of a complicat- 
ed equation by iteration. This may continue, 

183 D ' Blast JEtesponse of Airerafts Shulman s 130 hours 5 Porter » 2 hours 5 

VW1 B 45 minutes 

Two new programs (III, IV) were started. Program III is a modification 
of the previous programs to include a different forcing function* and to provide 
a means to stop the program when the beam resistance is zero 9 i 9 e» at the breaking 
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point o Trial runs xere started on this program , which is still in £ debugging 
stage o 

The other program started (IV) necessitated the use of the auxiliary drum 
for high-speed storage The forcing function for Program IV is the represents 
ation of the actual conditions and is expressed in terms of the Wagner function 
unlike the multi-linear approximations that were used in the previous programs 
(I II III). This required the use of the exponential library subroutine. 

Difficulty was experienced with the programming for the auxiliary drum 
but was overcome quickly . It is hoped that production runs can be started 
shortly . 

184 Do Scattering of Electrons frpm_Hy_droggn8 Newsteing .Iff, 464 minutes 

Final runs have been made on a series of programs. The results check 
to the required accuracy with the exact values which can be determined analytically 
for special values of the parameters. 

186 Co ■ Trackinjg_ResT3onse Characteristics of _jtheJ^gmaiL0Bg rator8 .Elkindg 

WWI„- 31 minutes 

During the past period yi minutes of computer time have been used. 
Most of this has been devoted to correcting errors made in previously computed 
power density spectra . Human operator transfer functions are in the process 
of being computed from the power density and cross pox^er density spectra 
determined by ¥$!. 

188 Co Effect of Gr avitx ...on Relator o '%ter Production in Oil Reservoirs 
Kerng Porter 8 2 hours i 'vl'jl, 33 minutes 

This problem is concerned with the solution of Laplace 8 s equation 
In two regions separated by a free boundary. Attempts to improve on an initial 
guess of the boundary position from the solution within the regions have not 
been successful . It is not known wheter this failure is due to inaccuracies 
in the solution (occasioned 'iyj an extremely low rate of convergence) or to a 
basic unsoundness of the method. An effort to increase the rate of convergence 
is being consides 



190 D. Zeeman and Stark effect in Positroniums Kendalls M£E„ 131 minutes 

The routine for computing the perturbed energy eigenvalues and partial 
lifetimes of the n' s 2 states of positronium has been completed and checked. 
.The program of computation is now underway. 

191 Bo -Earthquake Epicenter location by Gei ger^s M ethods Grine 9 .40 hours 5 

¥l£E 8 184 minutes 

The location of an earthquake epicenter observed by eight stations 
has been successfully determined. A generalisation of this program is now 
being developed so that data from any number of observing stations from three 
to forty may be fed 'in on a separate parameter tape. 
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195 Co Intestinal Motility; Farrar; Hamilton i WI 9 19 minutes 

Existing methods of analysing human and animal intestinal motility 
records do not permit an objective , accurate measurement of the total amount 
of motor activity presents, nor do they permit analysis of the component 
functions which make up the complex motility record. 

Intestinal motility in man is composed of several different basic 
wave forms , These includes (l) periodic n segmental" waves having a fre- 
quency of 8 - 11 e„p.m. (2) large "peristaltic" waves of low frequency 
and questionable periodicity 9 (3) random wave forms. Quantitative analysis 
of these components in a motility record has been impossible , as has been 
the evaluation of their relative contributions to the total motility record. 

During the past year a cooperative program has been undertaken be- 
tween the Gastroenterological section of the Evans Memorial Hospital and 
Mass. Institute of Technology. The aim has been to apply analytical methods 
used in the study of sound to the Interpretation of intestinal motility 
records. The tracings of rabbit ir.tra-enteric pressures and one- very 
simple human balloon-Kymographic record have -been analysed by generation of 
the auto-correlation function. -rain analysts has yielded numerical . values 
for the frequency and amplitude of the- varices motility components, as well 
as their relative contribution to the total record. The auto-correlations 
were initially obtained by using t -..;■■ mechanical correlator in Bldg. 32 of 
the Servome onanisms Laboratory, U /■■'.:■ over 9 this machine proved unsatisfactory 
for the analysis of the more ccmpl ..•:■■: human records since the wave excursions 
of the original data exceeded the n-.:-row limits of the machine. 

Generation of the autocorrelation function of a complex human 
motility record has been very satisfactorily achieved by the use of the 
Ifnirlwind I digital computer. Recordings are -made originally on rolls of 
paper having a 1 mm grid. Numerical amplitudes are then read at each of 
the equally-spaced 0- nua) points on the motility record. This data is pre- 
pared by Whirlwind personnel and the auto-correlation function calculated. 

It is thought that these exploratory studies done on Ttfhirlwind may 
form the basis for a more extensive study of intestinal motility by these 
analytical techniques . 

Dr. Farrar 

The work being done under problem .195 involves about 80 sets of auto- 
correlation data and Fourier transforms of the results. .The program being- 
used is the one developed by B In Ross under problem 107. 

Tape units have been assembled in groups of five data tapes with logging 
titles so that greater speed and a logging record can be. obtained. Three of 
these assembled tapes have been run. These indicated minor errors in the 
output routine o These errors are being corrected. Mass production of runs 

is expected soon. 

D „ . Hamilton 
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196. Single Address Computers Frankovich; 41 hours 5 Siegel, 36 hours ; 
Watson D 12 hours 1 l-MZ, 1429 minutes 

All modifications to the 1953 Summer Session computer (SS) to produce 
the 1954 Single Address Computer (SAC) wer > completed on schedule in time for 
use in the 1954 Summer Session Gcurs-? 6 531 « ^.n addition,, several existing 
errors were found in SS and corrected in S.\G„ The SAC computer has fixed 
point integer numbers only, but it has icsr 990 word pseudo-magnetic tape 
units and can combine numbers and instructions in ways forbidden in SS. A 
pseudo-time count is also made during SAC operations to inform programmers 
of the time a physically existent SAC would take to run their programs. 

2-97 . .T.bj"QQ_Mdr^ss_C_omputerg Demurjian s ? hours; Helwig» 40 hours; 
Bests TMI, 1530 minutes 

During the past several biweekly periods a hypothetical three-address 
drum computer (called TAG) has been developed. This was used by students dur- 
ing the earlier part of the 1954 summer session course 6.531« A brief descript- 
ion of TAC follows . 

TAG is a 3 address computer wiv,h 110 registers, of storage. Of these 9 
100 are contained on the drum (numbw-o-c; 00 Si .„»99) and 10 form a special high 
speed (zero access) storage (labelea :zQ 9 *. a9 z9) . 

Ik word in TAC is a series of 9 digits 8 letters or other characters 
(e.g., space „ tab, $ 9 ^ 9 etc.). R., ;: ,ter xO always contains 000000000 and 
cannot be altered by any instruction. Registers xl and x2 can automatically 
be dealt with as a double length numoer by many instructions. This double 
length location is labeled xx, 

TAG has four magnetic tape units « Each tape unit contains 40 blocks of 
information where each contains ten iAC words . 

Drum and tape speeds and operation times have been assumed for TAC. A 
hypothetical time, telling how long a program ran, is printed out each time a 
halt instruction occurs. Minimum l,.m,aney programming is possible. 

TAG words are typed in the form 
mn/abcdefghi 
and TAG tapes are terminated by 
ran/ start 

A conversion post mortem program detects ail the typing errors on TAC tape 
and prints' out the incorrect word (with the error in red) , the line on the 
print at which this oecured s and the distance along the paper tape. 

An operational post mortem prints' out the following sort of information 

when an alarm oc curs % e.g. 

stopped at 56l a 01020304 f 

0X5^99999999 02 I+99999999 03 j ^00000000 04J sq4q5060? 

Registers which have been changed s 
10 I ai01U213 ' 20 ) 4-23471211 
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Jump Table: 
(l,2,3,b) 2 t (6,7)3, 0,1 

The jump table expresses the path of control which the program followed . 
Loops are denoted by parentheses and exponents. Arbitrary loops within loops 
within , etc. are printed out by the program. 

198. Student Problems C oded for SAC ar n TAC : Staff; WI, 1175 minutes 

This problem number was set up for the us© of students enrolled in 
the special summer session course 6.531 (Digital Computers: Business Applicat- 
ions). 

Ml of the students in the course coded and debugged a problem assigned 
for solution on the simulated -three-address computer called TAC (described 
unpler problem 197) . ^He problem consisted of determining and printing out the 
day of the week corresponding to an arbitrarily specified date. 

Most of the students ; also worked on a second problem to be solved on the 
simulated single address computer SAC (see problem 196). For this second 
problem, the students were permitted to select" from among three suggested 
problems covering the maintenance of bank balance records, payroll calculations, 
and an autocorrelation calculation to' detect a seasonal sale variation «, respectively, 

199 • Laminar Boundar y 'layer of a Steady, Com pressible Flow in the .Entrance 
Region of a Tube s T-.oong; Wl, 97 minutes 

In connection with the research on heat transfer coefficients, 
recovery factors* and friction coefficients for supersonic flow of air in a 
tube, a theoretical investigation of -che characteristics of the laminar boundary 
layer in the entrance region has bean carried out. Partial differential equat- 
ions of continuity, momentum, and energy were developed, for the boundary layer. 
These were then transformed into a series of ordinary differential equations, 
to be solved by ¥KFX for several entrance Hach numbers and and thermal conditions 
at the tube wall. 

First, solutions -.re to be obtained for the case of constant viscosity 
and thermal conductivity. Then, the effects of temperature -dependence of these 
properties are to be tidied. 

Solutions of the first set of differential equations have been obtained 
for five cases of different entrance Mach numbers and thermal conditions at 
the tube wall. This represents approximately 5$ of the, entire job to be done 
by WHI. Routines are being prepared to solve the second set of differential 
equations. . 

The coding has been carried out in the K algebraic system" developed by 
Dr. H. Laning under S&EC problem 108. This system makes use of .Gill *s modified 
fourth-order Runge-Kutta method. 

200. .4.. Study of Recurrent Events s ■ • -3 . Jensen, ^0 hours; G:.;D-inaeen § 5. hours; 
Porter, 2 hours; TrftJI, 18 minutes ... 

An alternative method of discriminating ''between Bernoulli! sequences- of 
type S and Bernoulli! sequences of type N is by means 'of the density within 
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a given interval. If the interval is of length H 9 then all combinations of 
j9 Bernoulli! trials such that the number of ones is less than a threshold 
k compose our event, "Whenever recurrence time exceeds one 9 we say we detect 
a Bernouilli sequence of type S In terms of' density measurements, a sequence 
of type S is detected when the number of ones exceeds the threshold k. 
The program for this investigation has been coded and results obtained. 
Further investigations will be made to determine the effect of changing & ■ 
and k. 

201 o Study of the Ammonia Molecule s Meckler; I-JWI, 137 minutes 

Self-consistent wave-functions have been calculated for two geometric 
configurations of the ammonia molecule,, The integrals necessary for the final 
routine were derived by means of congruent transformations performed by .WI. 
These integrals were fed into the routine whtfch achieved self -consistent one- 
electron wave-functions in the framework of linear combinations of atomic 
orbitalso The results will be tabulated and discussed in the Quarterly Progress 
Report of the Solid State and Molecular Theory Group. 

205 o Electron Lattice Interactio n,. in. ..Solid as Schultz; W3X, 290 minutes 

To gain insight into electron- Lattice interactions in crystal when 
the coupling is too strong to be treated by perturbation methods (e.g. al- 
kali halides) exact results are beinr sought for a simple model consisting 
of two degenerate lattice oscillators 'corresponding to left and right going 
plane waves of the same wave length) Interacting with an otherwise fre© 
electron (moving in the direction el' "che plane xaves). 



matrices of order (n + l)(n+ 2)/2 arising .on neglect of all states 
having more than n phonons excited are geing diagonalized for n s l, 2,..-,. 6 
for divers values of coupling 9 oscillator wave lengths , and total momentum 
(electron plus wave field) with the following aimss 

1. to see row the effect of the interaction on the ground state depends 
on wave length 

2 to see if for sufficiently strong coupling, the lowest energy state 
has non-zero total momentum 

■ 3« "to provide- an accurate yarereiek for evaluating adiabatic methods used 
on real crystals and 

4. to provide a starting point for considering the effects of a periodic 
potential o 

Two values of the coupling constant (intermediate and strong) and two wave 
•lengths (intermediate and long) are being' investigated. For intermediate 
coupling good convergence has. been obtained. For strong coupling it is found 
that a better criterion than total number of quanta must be used in cutting 
off the infinite secular equation if good convergence is to be obtained. : Dr. 
Meckler is considering modifications in the present program which would include 
only those states with the n lowest diagonal matrix- components. 

206 . S3 r ectron_Ener??ies of the Molecule EZs Dalgarnoj TrJEE , 143 minutes 

The programs for the h varieties of coulomb integrals and the 3 varieties 
of . hybrid integrals were completed and tested. These programs have been used 
to obtain results for pertinent values of the parameters 9 thereby completing 
the problem (with the exception of the exchange integrals which will be done 
elsewhere) o 
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207. Check for RKACs Larson, 30 hours; Whalen, 30 hours; Porter, 3 
hours 5 WWI, 26 minutes 

The problem involves actually two check solutions for analogue machines 
for two systems of differential equations representing two different aircraft 
target configurations with control systems. The first is a" "16th order and the 
second is a 19th order system. Both are being solved by the 4th order Runge- 
Kutta method. The choice of interval is determined intuitively from past 
experience. Both routines have been punched on tape but no positive results 
have been obtained yet. 
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The following programs used computer time but did not reports 

An Interpretive Program 
Self-Consistent Molecular Orb rials 
Energy Bands In Crystals 
Synoptic Climatology 

Response of Mass Plastic Spring Sy 



Loading 

Utilizing a General Pun 
Switching-Circuit Desigi 



Digita 



suam to Transient 
1 Computer in 



Co Improved Power Spectrum Estimates 

B. Overlap Integrals of Molecular and Crystal Physics 

C. Tight Binding Calculations in Crystals 

D. Low Aspect Ratio Flutter 

B. Crosscorrelation of Blast Furnace Input-Output Data 

G. Eigenvalue • Problem for propagation of E„ M. Waves 

B. An 'Augmented Plane ¥ave Method As Applied "to Sodium 
Calculation of Vortical Antenna : Coverage Skeleton 



15 minutes 


36 minutes 


87 minutes 


^95 minutes 


8 minutes 


10 minutes 


4 minutes 


6 minutes 


. 13 minutes 


5 minutes 


16 minutes 


28 minutes 


419 minutes 


29 minutes 



lo3 Ojoeiating Sta tistics 

1„31 Cca puter_Time 

The following indicates the distribution of H£L time allocated to 
the S&EG Group during the four week period covered by this report* 

Programs 180 hours. 

Conversion 1 hour a 



Magnetic Drum Test 
Magnetic Tape Test 
Scope Calibration 
D emonstra t ions 
Total Time Used 
Total Time Assigned 
Usable Time B Percentage 
Number of Programs 



1 hour, 
1 hour;, 
1 hour 9 




o. mxmr&es 

13 minutes 
30 minutes 
.13 minutes 
16 minutes 



185 hours 20 minutes 
192 'hours, 57 minutes 
96.1$ 
722 
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•COMPUTER ENGINEERING 



• 2ol I'H'JI -System Operation 

(A. J. Roberts j L. L. Holmes) 

Computer dependability has been excellent. A failure of a relay- 
in the direct print out-punch out system resulted in the loss of 1 hour. 

.Computer shutdowns have been planned for the weekends of 11 and 18 
September. The first shutdown is to be used to realign our 600-amp alternator. 
The second shutdown will permit the relocation of the delayed printout equip- 
ment. Presently it is in P Row and is difficult to maintain. It will be moved 
to E.Row. 

(D, A. Morrison) 



The consolidated test. program 9 T- 3^2-16, is now in use in the ME 
marginal-checking routine. PMC can now be modified to' start over from a line 
number inserted in FF 4. A memo describing T-3^-32-16 is being written. 

3° ADMINISTRATION AM) P ERSONNEL 

Mew Non-Staff Personnel (R. A. Osborne) 

Sumner Bartlett is a new technician with the Construction Shop.' 
Burton Ewalt is a new member of the Drafting Department. 
Thomas J. MalIoy 9 Jr., is another new technician in the Construction Shop. 
Judith P. h/hite is a new Barta Building receptionist. 
Alfred W. Mojcicki has also joined the Construction Shop as a technician. 

Terminated Non-Stiff (PL A. Osborne) 

Diana Bierer 

4 a RECENT LIBRARY. ACQUISITIONS 

The following documents have bc;en recently acquired" "by the Barta Building 
Libra ry„ Room 1G9, 

No. Title Scutcb 



B 52 German-English Dictionary for Electronics Regen zxxL Regen 
Engineers and Physicists 

B 120 The Calculus of Finite Differences " L. IU halne-Theinson 

B 176 Petroleum Reservoir Efficiency and Well Committees of Standard 
Spacing . Oil and Humble' Oil 

B 19*r Kettenbruchen (Continuous Fractions) Oskar Perron 



DCL-13 



Page Ik 



C"152 Report to the Association for Computing Machinery 
First Glossary of Programming Terminology 

C-159 A Comparison of Machine Methods for Evaluating 
Certain Mathematical Functions 

C-l60 Proceedings? Railway Systems and Procedures Assn. 
. Inventory Management and Data Processing 

C-l6l A New Approach to Office Mechanizations Integrated 
Data Processing Through Common Language Machines 

C-162 Light on the Future 

G-I63 Electronics c o „ New Horizon in Retailing 

G-164 Integrating the Office for Electronics 

C-165 Short Range and Extended Forecasting by Statistical 
Methods; Air Heather Service Teeh. W Rept. 105-3? 

C-166 Colloquium on Linear Equations? Toeh. Kept. 35-5^ 

C-167 Use of Flight Simulators in the Design of Aircraft 
Control Systems 

C-168 - French York in the field of Digii&l Computation 

C-169 The York of the Swedish Board for Computing Machinery 

C-170 A High Speed Magnetic Drum Storage • 

G-171 WI Operation Logic 

F-185 Selective Analysis of the Scale of Atmospheric 
Turbulence (M.S. Thesis) 



A. CM. 
Y. Barkley Fritz 
Meeting— Chicago 
U.S. Steel Program 

I..B.M. 

Harvard Bus. School 
Am. Management Assn. 
Air Yeather Service 

O.N.R.L. 

H. Elmore Blanton 

ONBL-133-53 

OML Tech. Refit. 
OHR Tech. Memo. 
H-221 

J. .Howcroft-, 

J. "Smith 



